Localization of collisionally inhomogeneous condensates in a bichromatic optical lattice 
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By direct numerical simulation and variational solution of the Gross-Pitaevskii equation, we stud- 
ied the stationary and dynamic characteristics of a cigar-shaped, localized, collisionally inhomoge- 
neous Bose-Einstein condensate trapped in a one-dimensional bichromatic quasi-periodic optical- 
lattice potential, as used in a recent experiment on the localization of a Bose-Einstein condensate 
[Roati et al, Nature (London) 453, 895 (2008)]. The effective potential characterizing the spatially 
modulated nonlinearity is obtained. It is found that the collisional inhomogeneity has influence not 
only on the central region but also on the tail of the Bose-Einstein condensate. The influence depends 
on the sign and value of the spatially modulated nonlinearity coefficient. We also demonstrate the 
stability of the stationary localized state by performing a standard linear stability analysis. Where 
possible, the numerical results are shown to be in good agreement with the variational results. 

PACS numbers: 03.75.Lm,67.85.Hj,71.23.An 



I. INTRODUCTION 

Since Anderson predicted a localization of nonintcr- 
acting electron wave in solids with a disorder potential 
fifty years ago the Anderson localization has been 
observed and studied extensively in optics [H and acous- 
tics 3 and in Bose-Einstein condensates (BEC). In the 
study of Anderson localization in a BEC, disorder laser 
speckles 3 and quasi-periodic optical lattices (OL) [H 
have been used. Random speckle potentials are produced 
when light is refiected by a rough surface or transmit- 
ted by a diffusive medium Q. Billy et al. Q observed 
the exponential tail of the spatial density distribution 
when a ^^Rb BEC was released into a one-dimensional 
(ID) waveguide in the presence of a controlled disorder 
created by a weak laser speckle. A bichromatic OL is 
realized by a primary lattice perturbed by a weak sec- 
ondary lattice with incommensurate wavelength Q, and 
this system corresponds to an experimental realization 
of the Harper Q or Aubry- Andre model . Roati et al. 
@ studied the localization of a noninteracting ^^K BEC 
in a bichromatic OL. There have been many theoreti- 
cal studies of localization using the numerical solution of 
the Gross-Pitaevskii (CP) eq uation [l0| as well as using 
the Bosc-Hubbard model [Tl| in addition to the experi- 
mental studies under different conditions on disorder [T^ . 
There have been studies of localization in two- and three- 
dimensions and of the destruction of localization with 
the increase of nonlinear repulsion [13, [3] ■ 

A Feshbach resonance (FR) driven by a magnetic [l^ 
or optical field allows one to vary the atomic inter- 
action of a BEC in a controlled fashion [l^, thus cre- 
ating a noninteracting as well as a weakly interacting 
BEC for the study of Anderson localization. Fedichcv et 
al. [l^ predicted that the spatial variation of the laser 



field intensity by proper choice of the resonance detuning 
can lead to a spatial dependence of the atomic interac- 
tion, creating the so-called "collisionally inhomogeneous" 
BECs. The theoretical prediction has been demonstrated 
experimentally by Theis et al. with the ^'^Rb BEC [13 ■ 
Sakaguchi and Malomed studied the formation of solitons 
in such BECs [1^. There have been studies of matter- 
wave bright and dark solitons of the cubic-quintic nonlin- 
ear Schrodinger equation with time- and space-dependent 
nonlincarities [2lj in a collisionally inhomogeneous envi- 
ronment, and of dynamical effects of a bright soliton BEC 
with local and smooth space variations of the two-body 
atomic scattering length |22|. There have been studies 
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about how to introduce a space dependence in the non- 
linear interaction of the BEC in a controlled way 1231 . 
There have also been studies of soliton oscillations [2^ 
and dynamical trapping and transport [25| in collision- 
ally inhomogeneous BECs. 

Here we combine the two interesting settings, namely, 
the bichromatic OL and a collisionally inhomogeneous 
BEC (with a spatially modulated nonlinearity) , to study 
the statics and dynamics of a localized BEC in this set 
up. We assume that the spatial dependence of the non- 
linearity, induced by the external magnetic field of the 
OL, has the same form as the bichromatic OL, i.e., the 
nonlinear coefficient in the GP equation is proportional 
to the OL intensity [13, H^- We study the effects of 
the spatially modulated nonlinearity on the shape of the 
density envelope, and the stability of the stationary lo- 
calized states. The numerical results are shown to be in 
good agreement with the variational results, where appli- 
cable. On the other hand, the tail region of the stationary 
localized states is examined, where we find exponential 
decay in space indicating localization in a weak disorder 
potential. We study the location oscillation (oscillation 
of the center) and breathing oscillation of the localized 
states, induced by appropriate perturbation, employing 
numerical simulation of the GP equation. 

In Sec. |TT] we present a brief account of the ID GP 
equation and the bichromatic OL potential used in our 
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study, and a time-dependent variational analysis of the 
GP equation under appropriate conditions. We obtain a 
set of coupled evolution equations for the parameters of 
the localized state. The effective potential characterizing 
the spatially modulated nonlinearity is also described. In 
Sec. mil we investigate the effects of the spatially mod- 
ulated nonlinearity on the central and tail regions of a 
stationary localized BEG by a numerical solution of the 
GP equation using the split-step Fourier spectral method. 
For the central region, we demonstrate the stability of 
the localized states by performing a standard linear sta- 
bility analysis. In Sec. IIVI the oscillation dynamics of a 
collisionally inhomogeneous BEG in a bichromatic OL is 
studied numerically. In Sec. |V]we present a brief discus- 
sion and concluding remarks. 



II. ANALYTICAL CONSIDERATION OF 
LOCALIZATION 

We consider a cigar-shaped quasi ID BEG with inho- 
mogeneous atomic interaction described by wave func- 
tion u(x, t) satisfying the following dimensionless ID GP 
equation [27l. [28j 



du 



i CI u 



■ g{x)\u^u -I- V{x)u^ 



(1) 



with normalization jupda; = 1. The space variable 
X, time and energy are expressed in transverse har- 
monic oscillator units aj^ ~ ^/K/JrmJ), uj~^ and hw, 
where m is the mass of an atom and lo is the angular 
frequency of the transverse trap. As in the experiment 
of Roati et al. the potential V{x) is taken to be a 
bichromatic OL of incommensurate wave lengths: 



(2) 



with Al = 27r^.S(/A^, [l = 1,2), where A;'s are the wave- 
lengths of the OL potentials, sj's are their intensities, and 
ki = 27r/A; the corresponding wave numbers. In this in- 
vestigation, we take the incommensurate ratio of the two 
components ^ k2/ki = (^/5-l)/2 « 0.618033989.... In 
the actual experiment of Roati et al. however, the pa- 
rameter is fe/^i = 1.1972.... Without losing generality, 
we further take Ai = 10, and si = 10, S2 = 0.3si which 
is roughly the same ratio 82/ si as in the experiment of 
Roati et al. @. 

By means of the FR technique controlled by pro perly 
designed configurations of external optical fields [20, l2a] , 
the spatial variation of laser field intensity I{x) produces 
the spatial variation of the atomic scattering length. As 
the potential V{x) is also proportional to laser field inten- 
sity, it is reasonable to assume the spatial-dependence of 
the atomic scattering length is similar to V{x). Then, the 
atomic scattering length can be given as = aso+cV{x); 
here V{x) is the same as Eq. ([2]), a<,o is scattering length 



of the corresponding collisionally homogeneous system, 
and c is a constant coefficient related to the optical in- 
tensity and may be either positive or negative. Thus, 
the nonlinear coefficient g{x) in Eq. ([1]) has a spatial 
dependence of the form. 



g{x) = eo + e^Ai sm'^{kix). 



(3) 



The nonlinearity Eq is given by |27| eo = '2asoN/a± with 
the number of atoms, and e = 2cN/a± is the spa- 
tially modulated nonlinearity coefficient. Because the 
intensity-independent nonlinear coefficient gp may be al- 
tered independently of the other parameters |18| , in order 
to focus our attention on the effect of the spatially mod- 
ulated nonlinearity on the localization of the BEG, we let 
go = in the following. 

Usually, the BEG localized states formed on a bichro- 
matic lattice may occupy many sites of OL potential 
[13]. For certain values of the parameters, however, 
potential (O leads to localized states confined practically 
to a single site of the OL potential. When this happens, 
a variational approach for solving the GP equation is use- 
ful. To apply the variational approach to the localized 
BEG, we adopt the following variational ansatz 



u{x, t) 



ri/4 



exp 



2w^ 



X exp {i [7(x - xo) + I3{x - x^f (j)] ] ,(4) 



with w the spatial width of the localized state centered 
at a;o, 7 the linear phase coefficient, /3 the chirp and 4> the 
phase. These are time-dependent parameters [30| . The 
Lagrangian of the system is [3l|, |32| 



L{t) 



^ du du* 



dt 



dt 





du 




dx 



--9{x)\u\'-V{x)\u\' 



+N^Lm + NLi 



dx 
N / 1 



-f 7^ -f- 2p^w^ 



2 \2w^ 



(5) 



where the overhead dot denotes time derivative, the star 
denotes complex conjugation and 



L 



M 



4v27rw 



cos(2fcia;o) exp 



kfw^ 



1 ^ 

- ^ [1 - cos{2kixo) exp {-kfw'^)] 



1=1 



We use the variational Euler-Lagrange equation 

dL d dL 
da dt da 



0, 



(6) 
(7) 

(8) 
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where the variational parameters are cr = 0, xq, 7, /3, and 
w. The first variational equation using a = (j) yields 
N =constant. Wc take this constant to be unity and use 
it in the subsequent equations. The other choices of a 
respectively, lead to the following equations 



2V'2ttw 



Aiki sin(2fcia;o) exp 



1=1 



1 



Xq 

w 



Aih sin(2fciXo) exp (-/c, 

1=1 

■ 1, 

--2f3w = F{w,f3), 



2 2\ 



2^4 ^ 



1 dL 



M 



1 dLx 



w dw w dw 
The Hamiltonian of the BEC is 
dL -dL 



(9) 

(10) 
(11) 

= G(u;,/3).(12) 



H = xo 



dp 



\{^2+^" + '^P"^"]-LM^Ly. (13) 



Equations ([TT|) and determine the evolution of the 
width w once xo is known. Equations ^ and (|10p de- 
termine the evolution of center xq once the width w is 
known. To study the dynamics of the localized state, 
we insert Eq. © into Eq. ([TU)) to get an anharmonic 
effective potential Vcs'- 



cPxq 

VcSM 



■ = — T. \ l^cffM + Vcm ) , 

OXq OXq 

2 / ^ 

Ai cos(2fc;.To) exp ( -- 



4v27ni; 



(=1 



1 ' 

Kffv = -■^X!^'*^°^(^^'^o)'^^P(~^'^ 



(14) 

(15) 
(16) 



1=1 



The effective potential has two terms. The second term 
on the right of Eq. ([T4|) . 14fTv, arises from the bichro- 
matic OL and contributes to an attractive well, if \xo\ is 
small enough. The first term, VoffM, is induced by the 
spatial modulation of the nonlinearity and is actually a 
pseudo-potential (33l |. The pseudo potential is quasi pe- 
riodic and is a potential barrier or well depending on the 
sign and value of the coefficient e. 



III. STATIONARY LOCALIZED STATE 

The stationary states are obtained by setting the time 
derivative in Eqs. ([9]) — (jl2p to zero. Then the simplest 
solution of Eqs. ([9]) and (fTO]) is xq = and we consider 
below an immobile localized state at origin {xq = 0) with 
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FIG. 1: (Color online) The dimensionless effective potentials 
Veff "felt" by an immobile localized state for different e from 
Eqs. p5|) and p6p where the width w is calculated from Eq. 
HI. 



7 = 0. Equations (fTTj) and (|T2|) determine w and can be 
written as 



1 



0, 



(17) 



A, [1 - (1 + ^-^o^) exp f - 

2 

2«;4^A,fcfexp(-fcfu;2) ^0 



(18) 



1=1 



The effective potential felt by a stationary localized state 
at x = a;o = is obtained from Eqs. ([T4)) — p6)) and 
is plotted in Fig. [T] where the width w is obtained by 
solving numerically Eq. ([T8| . With the increase of £, the 
strength of both disorder and the quasi periodic effective 
potential increases. 

To understand the effects of the coefficient e, we obtain 
the stationary localized states by solving Eq. ([T|) numer- 
ically with real-time Split-step Fourier spectral method 
with a space step 0.04 and time step 0.0001. We checked 
the accuracy of the results by varying the space and time 
steps and the total number of space and time steps. To 
compare with numerics, the variational width is obtained 
by solving Eq. (HH). 

Typical numerical and variational results for density, 
width and Hamiltonian of the localized state at a; = 
are exhibited in Fig. [2] for — 1 < e < 1. The variational 
Hamiltonian is obtained from Eq. p3l) and the numer- 
ical Hamiltonian is obtained from H ~ J^^[\u'\'^ /2 + 

g{x)\u\^/2 + V{x)\u\'^]dx. Figures [2] (a) and (b) exhibit 
the density of the stationary localized states correspond- 
ing to £ = — 1 and 1, respectively. Figures [2] (c) and (d) 
exhibit the variation of w and H with e. The numerical 
width w in Fig. [5] (c) is \/2 times the root-mean-square 
(rms) size of the BEC. Figures [2] indicate that the numer- 
ical results are in good agreement with the variational 
results. Figure [2] (c) shows that the width decreases and 
Fig. [11(d) shows that the Hamiltonian increases with the 
change of s from negative to positive. The dependence 
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FIG. 2: (Color online) Numerical (lines) and variational 
(chain of symbols) densities |up of the localized BEC ver- 
sus X for (a) e — —1, (b) e = 1. (c) Numerical (solid lines) 
and variational (dashed lines) dimensionless widths w versus 
e and (d) Hamiltonian versus e. (All quantities are dimen- 
sionless.) 



of the variational width on e can be quahtativcly under- 
stood as follows. The height of the central well of the 
effective potential increases with increasing e as shown 
in Fig. [T] Hence the central part of the localized state 
with a Gaussian shape becomes narrow with the increase 
of e, as can be seen from Figs. [2] (a) and (b). The vari- 
ational Gaussian ansatz only represents this central part 
and hence the variational width decreases with increas- 
ing s. However, the numerical width (rms size) shown 
in Fig. [5] (c) receives nontrivial contributions from both 
the central Gausian part and the extended exponential 
tail of the wave function (viz. Fig. 2] (a)), making it 
difficult to predict even qualitatively the variation of the 
numerical width with e. It is interesting to compare these 
results with those of a collisionally homogeneous conden- 
sate where a constant negative (attractive) nonlinearity 
leads to a reduction of the width and a constant pos- 
itive (repulsive) nonlinearity increases the width of the 
localized state [l^. The numerical width is larger than 
the variational width in Fig. [5] (c) due to the long expo- 
nential tail of the numerically obtained BEC (viz., Fig. 
21). Also, the difference between the variational and nu- 
merical widths increases for larger values of e, because for 
larger e the exponential tail is more pronounced resulting 
in a larger width (viz.. Fig. 

It is important to investigate the stability of the sta- 
tionary state governed by Eqs. pT)) and ([T^ against 
perturbation by a standard linear stability analysis. In- 
troducing small fluctuations around the stationary solu- 
tion (wo,/3o), w'{t) = w{t) - wo,P'it) = Pit) - /3o, and 
linearizing Eqs. pT|) and ([T^ indicated by them, a set 



of two linear equations are obtained 
dw'{t) 



dt 
dp'{t) 
dt 



= F^{wo.Mw\t)+Fp{wo,fio)P\t), (19) 
^ G^(«;o,/3oKW + G^K,/3o)/3'W, (20) 



where the subscripts w and /3 denote a derivative with 
respect to the respective variable. Assuming the solution 
of w'it) and /?'(<) m exponential form, ^ exp(£<), the 
eigen value £ is 



2£ = F^{wo.Po) + Gp{wo,Po) 

±{[i^»K,/3o) +G;3K,/3o)] 

-l-4F0(u;o,/3o)G^(w;o,/3o)} 



1/2 



(21) 



From Eqs. ([TT]) and ^ , we find F^iwQ.ISo) = 
0, G/3(wo,/3o) = 0, and 

Fp[wo,(io) = 2wQ, (22) 
2 ^ 

Gw[wq,(3q) = g + 2u;o^^ifcf exp [-kfwl) 



1=1 



+ 



-3 + (3 + 2kfwl + kfwl) 



X exp 



- ; 2 2 



which leads to the eigen- values 

£ = ±[F0{wn,PQ)G^{wo,P^f'^. 



(23) 



(24) 



To investigate the stability, Eq. (|18p is first solved to get 
as a function of e. This result is then inserted in Eqs. 
(p^l) and ^ to get F/3(wo,/?o) and G^(wo,/3o). The 
graphical representation of the two functions is shown in 
Figs. [3] (a) and (b). In the case of a small coefficient e, 
we can find that Gu,(wo, Po) < and Fp{wo, Pq) = 2wo > 
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FIG. 4: (Color online) (a) The numerically obtained dimen- 
sionless density |u|^ versus dimensionless x for different e. (b) 
The dimensionless localization lengths L versus e. The local- 
ization length L is obtained by exponential fitting to the tails 
of density distribution with ~ exp( — |a;|/L) 



0. Thus, both the eigen-values from Eq. ([M)) must be 
imaginary, and the localized state from Eqs. (|17p and 
([TS]) is stable against small perturbation. 

Anderson localization in a weakly disordered potential 
is characterized by a long exponential tail of the local- 
ized state. For coUisionally homogeneous condensates, 
the experimental [sj and theoretical [l^ investigations 
have demonstrated that the localized BECs have an ex- 
ponential tail for wcakly-intcracting or non-interacting 
BEC in a quasi-periodic OL. In order to observe the ef- 
fect of spatially modulated nonlinearity on the tail region, 
we plot in Fig. S] (a) the density distribution |up of the 
stationary BEC on log scale. As we see, the long expo- 
nential tail extends from x « ±2 to x ~ ±20, whereas the 
central part of density for |a;| < 2 represents a Gaussian 
distribution. By an exponential fitting of the exponential 
function exp(— |a;|/X) to the tails of density distribu- 
tion, the localization length L versus e is illustrated in 
Fig. m (b) which shows that L increases nonlinearly with 
e. An increase in e represents a decrease in disorder thus 
resulting in larger values of localization length. 



IV. DYNAMICS OF LOCALIZED STATE 

To get further insight into the effects of the spatially 
modulated nonlinearity e on the localized states, we now 
study some dynamics of the localized state. First, we 
study numerically the oscillation of the localized states 
in an external potential. According to Eqs. pl| - 
the motion of the localized BEC can be approximately 
regarded as that of a particle inside an effective potential 
VcS- Because of the exponential tail and elasticity of the 
localized state, although the variational results may not 
be good for the dynamical evolution, they can provide a 
qualitative physical understanding of the dynamics using 
the effective potential. To study the motion of the local- 
ized state, first we create a stationary localized BEC with 
spatially modulated nonlinearity in the bichromatic OL. 
Successively, at i = 0, we suddenly introduce an initial 
momentum po = 0-1 by u{x) — >■ uo(a;) exp(ipoa;), where 
uq is the wave function of the stationary localized BEC. 
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FIG. 5: (Color online) The center of the localized state xo 
versus time t (both dimensionless) during the location oscil- 
lation initiated by suddenly introducing an initial momentum 
po = 0.1 by the transformation m — >■ uq &q>{ipox) for e — (a) 1, 
(b) 0, and (c) —1: numerical result (full line) and variational 
result (dotted line). 



From the experimental point of view, the initial momen- 
tum can be given by suddenly moving the OL 



3J. It is 



found that after the perturbation, the density envelope 
suffers an abrupt change but remains localized. Actually, 
the localized BEC is an elastic object and not a rigid one. 
Hence, both the center and the density distribution of the 
localized state perform oscillations. 

The evolution of the center xq of the localized BEC as 
obtained from numerical simulation (full line) is shown in 
Fig. [S]where xq is obtained by instant Gaussian-function 
fitting to the central region of the density distribution. 
The variational results are also shown in this figure. The 
top, middle and bottom panels correspond to e = 1,0, 
and —1, respectively. We find that, in general, the oscil- 
lation of the localized state could be quasi periodic after 
an initial damping, which can be explained on the ba- 
sis of energy conservation. Because of the deformation 
of the density envelope, a part of the initial kinetic en- 
ergy is converted into elastic energy of deformation and 
this leads to the damped oscillation. This deformation 
will be larger for large e(= 1), when the localized state 
is loosely bound with large exponential tail (viz. Fig. |4] 
(a)). Consequently, an oscillation of the center xq with 
rapidly varying amplitude is found denoting easy periodic 
transfer of energy between location oscillation and defor- 
mation. For small £(= —1), the localized state is more 
compact and tightly bound, so that it can be treated like 
a rigid object. The exchange of energy is less probable in 
that case, and a periodic oscillation of the center xq with 
constant amplitude is found. Eventually, the energy of 
deformation is liberated leading to an increase in the am- 
plitude of oscillation. During the subsequent oscillation 
cycle, the conversion between the kinetic energy and clas- 
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solved to obtain the oscillation of the central part and 
this result is shown in Fig. [HI (d) for e = 0. (The Gaus- 
sian variational ansatz without any exponential tail only 
carries information about the central part.) There are 
two regions of the localized state which oscilate with two 
distinct frequencies: (i) the central region with Gaussian 
distribution and (ii) the outer tail with exponential dis- 
tribution. The net result is the harmonic oscillation of 
the width with a small frequency and large amplitude due 
to the oscillation of the exponential tail, which is modu- 
lated by rapid oscillation coming from the central Gaus- 
sian part. The variational frequency as obtained from 
Fig. [6] (d) is found to be identical with the frequency of 
modulation of the numerical width, which confirms that 
the modulation in Figs. in](a), (b) and (c) is coming from 
the oscillation of the central part of the condensate. The 
renewed oscillation in Figs. [5] and |6] confirms the stability 
of the stationary localized BECs. 



FIG. 6: (Color online) Numerical (N) result for dimensionless 
pulse width w of the localized state versus dimensionless time 
t during breathing oscillation initiated by suddenly changing 
the strength of the secondary lattice S2 from 3 to 3.5 for e — 
(a) 1 , (b) , and (c) —1. We also show the variational result 
(V) of pulse width w from a solution of Eqs. and (|12|) 

with condition w{t = 0) = 0.9688 and l3(t = 0) = 0. 

tic strain energy causes the quasi-periodical movement of 
the localized state. In the variational formulation the ex- 
change of kinetic energy to the energy of deformation is 
not allowed and the resultant oscillation is of a fixed am- 
plitude without damping. Nevertheless, the numerical 
frequency of location oscillation is in agreement with the 
variational frequency within an estimated error of about 
2.5%. As pointed out in Sec. IIIIl a positive e leads to 
a tighter trapping and vice versa. Then, with the same 
initial velocity, the tighter trapping causes the localized 
state to oscillate with larger frequency and smaller ampli- 
tude, and a weaker trapping leads to a smaller frequency 
and larger amplitude of oscillation as indicated in the 
numerical results of Fig. [S] The variational frequency 
follows the same trend as e is changed from positive val- 
ues to negative values. 

Next wc consider a breathing oscillation of the local- 
ized BEC, started by suddenly changing the strength of 
the secondary lattice S2 from 3 to 3.5 at t = 0. We in- 
vestigate how the breathing oscillation of the localized 
BEC is changed by the spatially modulated nonlinear- 
ity. Then, the nonlinearity in Eq. ([T]), g{x), also changes 
with the new OL potential. We present numerical results 
in Fig. ini (a) , (b) , and (c) for the time evolution of the 
width w of the localized state for e = 1,0 and —1, re- 
spectively. The variational equations PT|) and were 



V. SUMMARY 

In this paper, using the numerical and variational solu- 
tion of the time-dependent GP equation, we studied the 
stationary and dynamic characteristics of a cigar-shaped 
localized BEC with spatially inhomogeneous nonlinearity 
in a bichromatic quasi-periodic ID OL potential. This in- 
vestigation reveals that the spatially inhomogeneous non- 
linearity produces a pseudo-potential which changes the 
strength of the disorder and the height of the quasi peri- 
odic effective potential felt by the localized BEC. With a 
larger spatially modulated coefficient e, the localization 
length and Hamiltonian will be larger. 

We also study the stability of the stationary localized 
state using the linear stability analysis and find it is dy- 
namically stable under small perturbations. The stabil- 
ity is also verified by numerical simulation. In respect 
to dynamics, wc investigate the location oscillation (os- 
cillation of the center) and breathing oscillation of the 
localized BEC, and find that both oscillations are quasi 
periodic because of the quasi periodic effective potential. 
The frequency of quasi periodic oscillations of the center 
of the BEC increases as e increases. For the breathing 
oscillations, the two exponential tails also are symmetric 
around the center at x = 0. The present study is use- 
ful for an understanding of the statics and dynamics of 
Anderson localization and for planing new experiments 
with coUisionally inhomogeneous BEC. 
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